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ABSTRACT 
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It is argued that the nanosecond giant subpulses detected recently in the Crab 
; pulsar are generated by means of the coherent curvature radiation of charged 

j> ; relativistic solitons associated with sparking discharges of the inner gap potential 

| Q ■ drop above the polar cap. 
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1. Introduction 



Although only four pulsars are known to emit giant pulses (Lundgren et a. 1995; Cog- 
nard et al. 1996; Romani & Johnston 2001; Johnston & Romani 2003), understanding the 
mechanism of their radiation can potentially lead to understanding the longstanding prob- 
lem of a pulsar radio emission. Recently, the detection of extremely short and powerful, 



2 ' 2-nanosecond - 1000 Jy, subpulses within the radio giant pulses from the Crab pulsar has 
been reported by Hankins et al. (2003, HKWE03 hereafter). This is an observational re- 
sult of extraordinary importance, for it can shed a new light onto the mystery od coherent 
pulsar radio emission. HKWE03 argued that these nanosecond giant subpulses were true 
temporal modulations associated with an explosive collapse of nonlinear plasma turbulence 
(Weatherall et al. 1998). In this letter we argue that the time-scale of the Crab pulsar giant 
nanosecond subpulses is consistent with the angular beaming due to the curvature of dipo- 
lar field lines (different from the conventional angular beaming due to the pulsar rotation), 
thus with the coherent curvature radiation. Furthermore, we demonstrate that the observed 
fluxes of these giant subpulses are consistent with the emission of charged relativistic soli- 
tons, generated by means of the modulational instability of the strong Langmuir turbulence 
associated with a sparking discharge of the pulsar's polar gap (Melikidze, Gil & Pataraya 
2000, MGP00 hereafter; Gil, Lyubarski & Melikidze 2004, GLM04 hereafter). 
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2. The time-scales 

2.1. The Alignment Time-scale 

Let us consider a number of localized (point-like) sources of the broadband coherent 
curvature radiation, moving relativistically (with the Lorentz factor 7 = (1 — /3 2 )" 1 / 2 , where 
(5 — v/c ~ 1) along a narrow bundle of dipolar magnetic field lines with a radius of curvature 

p = l.Sr^VHO 7 cm, (1) 

where s = d/r p is the normalized polar coordinate of the central line of the bundle, r 6 = 
r cm /R is the normalized emission altitude r (for the frequency v = 5 GHz, 5 < r 6 < 35 in 
the Crab pulsar, Kijak & Gil 1998), d is the distance from the dipolar axis to the foot of a 
dipolar field line, r p « 10 4 P -1 / 2 is the polar cap radius and R = 10 6 cm is the neutron star 
radius. The geometry is schematically illustrated in Fig. 1 where a narrow bundle of field 
lines is marked. The curvature radiation is emitted along dipolar field lines into a narrow 
cone with an opening angle 9 = l/j (e.g. Jackson 1975). In what follows we assume that the 
perpendicular dimension of each source measured by the angular extent of the bundle of field 
lines does not exceed 9 (validity of this assumption is justified later on in the paper). We also 
assume that other bundles carrying sources of coherent curvature radiation are separated by 
at least 9, thus only one narrow bundle contributes to the radiation observed at a given 
rotational phase = 2-irt/P, where t is the intrapulse time. For a fixed arrangement of the 
observer's direction, the field line polar coordinate s and the frequency dependent emission 
altitude r^u), each source is aligned within I/7 with the observer during the time interval 
5t = l/v = p9/v = p7 _1 /v (see Fig. 1), where p is described by equation (1). Therefore, a 
continuous stream of relativistic sources moving along a curved trajectory would illuminate 
the observer during the time interval 

At = P1 - l /c , (2) 

provided that this "alignment" time-scale is longer than the conventional angular beaming 
time scale At' = 7~ 1 P/27r related purely to the pulsar rotation. In the case of the Crab 
pulsar this means that p > 1.6 x 10 s cm in the case of the Crab pulsar. HKWE03 concluded 
that the conventional angular beaming cannot explain the 2-nanosecond duration of their 
giant subpulses, since it indeed requires 7 ~ 10 6 , much above the expected values 7 < 1000. 
Since, incidentally, in the case of the Crab pulsar At' ~ At, their conclusion applies also to 
the alignment time scale related to the field line curvature (eq. [2]). 
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Fig. 1. — Geometry of the observed curvature radiation emitted by the source (S) moving 
relativistically along a narrow bundle od dipolar field lines with the radius of curvature p 
(see text for explanation of other symbols). 
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2.2. The Apparent Time-Scale 

The situation can be quite different if the stream is not continuous but some sources are 
distinguished by emitting slightly more power V than the average. From the observational 
point of view it means that there are very few sources emitting towards the observer. For 
simplicity let us first consider a single source, which is marked in Fig. 1 in three different 
positions A, B and C. Following Jackson (1975) we can argue that the duration of the 
observed impulse related to the curvature of field lines is much shorter than that predicted 
by equation (2). Let us consider two particular points along a given dipolar field line: 
first (A) at which the source Si becomes aligned (ray 1) and second (B) at which it becomes 
misaligned (ray 2) with the observer O (Fig.l). The source covers a distance Z = p/7 between 
these points during a time interval 5t = P7 _1 /v. During this time the radiation emitted at 
the first alignment point (A) travels a distance L = c5t = (c/t>)p7 _1 = p/((3~f). Thus, the 
radiation overtakes the source only by a distance A = L — I = (1/(3 — l)(p/^) p/7 3 (when 
the radiation emitted at point A reaches point B, the source of this radiation arrives at point 
C) . This is the pulse length in space, and thus the duration of the observed impulse emitted 
by the considered source is 

Ar = ^ = 4 = 7 - 2 At, (3) 
c cy 

(Gil 1985) where At is the alignment time scale (eq. [2]). Therefore, if the longitudinal (along 
field lines) source dimension is smaller than A pa p/7 3 , then the apparent duration of the 
impulse At is 7 2 times shorter than the actual time of alignment of the source with the 
observer. It is important to emphasize that this apparent shortening effect has nothing to 
do with special relativity (time dilation) and/or Doppler effects, as both At and Ar are 
referred to the same observer's frame. This effect occurs for any localized source of emission 
(radiation, particles, etc.) moving towards the observer (target) with a speed slightly lower 
that the emitted waves (particles). For example, one can consider an acoustic analog in 
which a highly directional loud-speaker is moving along a curved trajectory with a velocity 
close to the speed of sound. The duration of the acoustic impulse is much shorter than the 
actual time of geometrical alignment of the beam pattern with the recording device. The 
acoustic analog of equation (3) is Ar = 2 At, where 6 is the beam-width of the emitted 
pattern (in radians). 

Both isolated 2-ns giant subpulses, as well as a sequence of a number of such subpulses 
(when unresolved this sequence constitutes a normal giant pulse in our view), can be seen 
in Fig. 1 of HKWE03. Within our model such case corresponds to a sequence of sources, 
each having a longitudinal dimension smaller than A and separation between the adjacent 
sources larger than A, where A ps p/7 3 is the spatial length of the impulse associated with 
each source. For each giant subpulse detected by HKWE03 the observed time duration 
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At = 2 x 10 9 s (eq.[3]), which leads to the condition 

Ps » 0.6 72 3 , (4) 
where p 8 = p/10 8 cm (eq.[l]) and 72 = 7/IO 2 . 

3. Energetics 

3.1. Brightness Temperature 

As reported by HKWE03, the ultra-short nano-giant subpulses detected at v — 5 GHz 
in the Crab pulsar often exceed fluxes S = 1000 Jy = 10~ 20 erg s _1 Hz _1 cm _1 . These fluxes 
can be converted into the brightness temperature Tj, = SD 2 /[2k(uW) 2 ], where D = 2 kpc = 
6 x 10 21 cm is the distance to the pulsar, k is the Boltzmann constant and W is the time scale 
of the corresponding emission process (e.g. McLaughlin and Cordes 2003). Assuming that 
W = 2 x 10 _9 s (HKWE03) one obtains the extraordinarily high brightness temperatures 
T b ~ 10 38 K, implying by far the most luminous emission from any astronomical object. 
However, within our model the actual emission process corresponding to a giant subpulse of 
duration At = 2 x 10~ 9 s occurs over a much longer "alignment" time interval W — At — 
■y 2 At (eqs. [2] and [3]), which leads to 7 4 times lower brightness temperatures. Since 7 is of 
the order of 100, then Tf, < 10 30 K, consistent with normal giant pulses in the Crab pulsar 
and other pulsars (see Fig. 1 in McLaughlin and Cordes, 2003). 



3.2. Luminosity and Power 

Let us then consider a source(s) moving along a dipolar field line (Fig. 1) and emitting 
a coherent curvature radiation with an intrinsic power V . While moving over the alignment 
distance I = p7 _1 , the emitted energy E = VAt, where the alignment time interval At = 
P7 _1 /c (eq.[2]). This energy is received by the observer in a much shorter time At = 7 _2 At 
(eq.[3]), and therefore the apparent luminosity 

£=^=7 2 P. (5) 

Thus, the nanosecond giant subpulses appear 7 2 = 7! 10 4 more luminous than the intrinsic 
power of their source(s). Also the apparent fluxes S oc C are 7 2 times overestimated. It is 
important to emphasize that this kinematical boosting is a direct consequence of an apparent 
shortening of the impulse duration described by equation (3) and again has nothing to do 
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with special relativity and/or Doppler effects (Ginzburg 1979). According to the discussion 
below equation (3), an acoustic analog of equation (5) is £ = 0~ 2 V. 

The apparent fluxes S ~ 1000 Jy of the nanosecond giant subpulses can be formally 
converted into the emitted luminosity C = S0 2 D 2 Au. With D = 2 kpc, 6 = 7 5 T 1 10~ 2 and 
Av ~ 10 9 Hz one obtains C ~ 72" 2 10 28 erg/s. This is an extremely large value, compara- 
ble with the total radio luminosity of the Crab pulsar Lr = 3.5 x \0 25+x erg s _1 , where 
x = logL ps 3.6 mJy kpc 2 from the Pulsar Catalog (Taylor, Manchester & Lyne 1993). The 
existence of such powerful localized sources of coherent radio emission seems highly question- 
able (of course, this argument is not independent from the extraordinarily high brightness 
temperatures T b ~ 10 38 K derived in Section 3.1). However, according to equation (5), the 
actual power of the emitted radiation V = r )~ 2 S9 2 D 2 Av = 7~ 2 £, which for the parameters 
given above yields approximately 

V ~ 7 2 ~ 4 1 24 erg s _1 . (6) 

This is still a very large power (consistent with T b < 10 30 K) but much smaller than the total 
radio luminosity of the Crab pulsar. In the next section we discuss a mechanism of coherence 
of the curvature radiation that can produce shots of radio emission of nanosecond duration 
with the intrinsic power corresponding to the values described by equation (6). 

Finally we can check whether our localized source can be supplied with enough kinetic 
energy to emit giant subpulses. To estimate the energy of the source let us note that the 
product of plasma number density and the cross-section of the corresponding flux tube is 
a constant value equal to Kn GJ h 2 , where n G j is the Goldreich & Julian (1969) density, k is 
the Sturrock (1971) multiplication factor and h is the gap height (equal to the spark size). 
Assuming that the perpendicular source dimensions are determined by the spark size, we 
can write that the kinetic energy associated with the source is E ST ps nGjKmc 3 h 2 A/c. On the 
other hand, the giant subpulse energy E GP = At = ^y 2 VA/c (eqs.[3] and [5]). Obviously 
E GP should be smaller than E sv , which leads to the condition > 500, consistent with an 
independent estimate given by GLM04. 

4. Coherent curvature radiation 

MGP00 proposed a model for generation of the coherent pulsar radio emission, based on 
the idea of the polar gap discharging via a number of localized sparks (Ruderman & Suther- 
land 1975). The sparking phenomenon creates a succession of plasma clouds moving along 
dipolar magnetic field lines. The overlapping of charged particles with different energies 
from the adjacent clouds ignites a strong Langmuir turbulence via the two-stream instability 
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(Usov 1987; Asseo & Melikidze 1998). In the pulsar magnetospheric conditions, this tur- 
bulence is subject to modulational instability, which leads to the formation of "bunch-like" 
charged solitons capable of generating the coherent curvature emission at radio- wavelengths 
A. Thus, the condition A < A = p/7 3 = (p 8 /7|)10 2 is naturally satisfied at centimeter 
wavelengths (note that according to eq.[4] A = 60 cm, which of course corresponds to 2 
light nanoseconds). Additionally, the angular extent Aijj of the bundle of dipolar field lines 
associated with a given spark should not exceed 9 = 7^. It is easy to show that in the Crab 
pulsar Aip ~ 0.001 < 7- 1 (for details see Gil k Sendyk 2000). 

MGP00 calculated the power of the soliton curvature radiation in the vacuum approx- 
imation, that is without taking into account the influence of the magnetospheric plasma. 
GLM04 reconsidered the curvature radiation of a point-like charge moving through electron- 
positron plasma. They demonstrated that the radiation power is suppressed by a factor of 
about 10~ 2 - 10~ 3 but still at a high enough level to explain the observed pulsar luminoc- 
ities. Applying the results of MGP00 modified by GLM04 to the case of the Crab pulsar, 
we obtain that the power radiated by one soliton can be as high as L\ ~ 10 21 erg/s. Thus, 
for an incoherent superposition of N sources of the coherent radio-emission (solitons), each 
smaller than A and separated from each other by more than A, we have 

L ~ A^IO 21 erg/s, (7) 

(for details see eqs. [14] - [18] in MPG00). Comparing equation (7) with the power required 
for the Crab giant subpulses expressed by equation (6) one obtains the Lorentz factors 
72> lO^/iV 1 / 4 , which is realistically about 3-4 (consistent with an independent estimate for 
7 ~ 300 — 400 given by GLM04). Inserting 72 = 3 into equation (4) we obtain p 8 > 16, 
which, according to equation (1), means r^s" 1 ~ 90. Since r 6 < 35 (Kijak & Gil 1998) 
this implies that s < 0.07. Remembering that s is a normalized polar coordinate (0 - for 
the pole and 1 for the polar cap edge), we conclude that only those field lines that originate 
very close to the magnetic axis can be involved with the nanosecond giant sub-pulses. Such 
field lines can be aligned with the observer in the radio emission region only within a very 
narrow range of longitudes near the so-called fiducial plane, containing both the spin and 
magnetic axes. It is natural to associate this plane with the peak of the pulse profile. Thus, 
within our model, the giant subpulses should occur only within a very narrow range of 
longitudes near the peak of the radio profile. To the best of our knowledge, giant subpulses 
in the Crab pulsar seem to occupy the narrow range of phases within about 1% of the pulse 
window near the peak of the main-pulse (HKWE03). This is consistent with the phase range 
2^/360° = s(2.°4/360°)(r 6 /P) 1 / 2 , which for s < 0.07 and r 6 < 35 (see above) is smaller than 
0.015. 
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5. Discussion 

In this letter we argue that the 2-ns giant subpulses detected in the Crab pulsar are 
due to the curvature radiation of one or at most several solitons associated with a sparking 
discharges occurring near the local surface magnetic pole. This is the most elementary emis- 
sion event that can be observationally resolved. The curvature radiation from the particular 
soliton(s) can be observationally distinguished by means of the kinematical boosting (eq. [3]) 
only if they are slightly more powerful than other solitons. This implies that the sparking 
event associated with giant subpulses is more energetic than an average one. Both luminos- 
ity and characteristic frequency of the soliton curvature radiation strongly depend on the 
Lorentz factor 7, which in turn depends on the accelerating potential drop above the polar 
cap. Interestingly, pulsars exhibiting giant pulses are distinguished not only by their bright- 
ness, but also by extremely high values of the so-called complexity parameter a = (r p /h), 
where r p is the polar cap radius and h is the polar gap height (Gil & Sendyk 2000). Since 
the accelerating potential drop Ay oc h 2 and h < r v j\[2 (Ruderman & Sutherland 1975), 
then pulsars with very high values of a » 1 must have a large reservoir of the maximum 
available potential drop over the actual potential drop (AV max / AV = a 2 /2). This reservoir 
can be occasionally used to create exceptionally energetic spark(s). 

The values of the complexity parameter a can well exceed 100 in normal pulsars, while 
in millisecond pulsars they are limited to about 40 (see Fig. 1 in Gil & Sendyk 2000). 
Therefore, when pulsars in both groups are ranked with respect to a, pulsars with giant 
pulses occupy the top of the lists (Table 1). Also, the Vela pulsar, which was reported to 
exhibit some kind of giant pulse behavior (Johnston et al. 2001) is very high on the list. 
Moreover, sporadic large amplitude pulses (LAP) from two millisecond pulsars (J1959+2048 
and J0218+4232) have been reported very recently by Joshi et al. (2003). It should be 
emphasized that pulsars showing giant pulses have also high values of the magnetic field 
B LC > 10 5 G at the light cylinder. This might, however, be a coincidence because roughly 
-Blc ~ a2 /Pi where P is the pulsar period (Table 1). Although B LC is a good parameter 
to make a list of giant pulse candidates, only a value of a 2 has a physical meaning within 
our model. In Table 1 we propose a number of candidates for giant pulses with high values 
of the complexity parameter (a > 10 for millisecond pulsars and a > 90 for normal pulsars) 
and relatively low values of -Blc- The detection (non-detection) of giant pulses from these 
candidates would confirm (refute) our scenario. PSR Jll 19-6127 with low -Blc ~ 5 x 10 3 G 
seems the most interesting case. 

The case of PSR J1959+2048 deserves some additional discussion in the light of the 
model discussed here. As emphasized by Johnston & Romani (2003), this pulsar is ex- 
ceptional in the sense that despite the high value of B LC it does not show giant emission 
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Table 1. Pulsar parameters 



PSR J 


P 


P 


a 


Blc 




(msec) 


(io- 15 ) 




(10 5 G) 


0534+2200 


34 


420 


249 


8.7 


(2) 


0540-6919 


51 


479 


198 


3.3 


(6) 


1513-5908 


150 


1540 


138 


0.4 


(21) 


1124-5916 


135 


745 


120 


0.4 


(22) 


1617-5055 


69 


137 


114 


0.8 


(12) 


1420-6048 


68 


83 


99 


0.7 


(14) 


1119-6127 


408 


4002 


95 


0.05 


(75) 


0835-4510 


89 


125 


94 


0.4 


(20) 


1824-2452 


3 


0.0016 


33 


7.3 


(4) 


1823-30A 


5.4 


0.0034 


28 


2.4 


(8) 


1939+2134 


1.6 


0.0001 


23 


8.8 


(1) 


0218+4232 


2.3 


0.00007 


16 


3.1 


(7) 


1959+2048 


1.6 


0.00002 


14 


3.6 


(5) 


2129+1210E 


4.6 


0.00018 


14 


8.2 


(3) 



Note. - Pulsars marked in boldface show gi- 

ant pulses, others are candidates for giant pulse emis- 
sion. Numbers in parenthesis denote ranking position 
in B LC = 2850P a5 P- 2 ' 5 . The complexity parameter 
a = 425(P) 2 / 7 P- 9 / 14 . 
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(although Joshi et al. (2003) found one significant LAP (out of million pulses) with an in- 
tensity 129 times the mean intensity). This seems to be consistent with our model since 
this pulsar (which is fourth ranked with respect to 5lc) has a relatively small value of the 
complexity parameter a as compared with two other millisecond pulsars showing giant pulses 
(Table 1). In fact, since AV^ax/AV = a 2 /2, the smaller the value a, the smaller the reservoir 
of the polar gap energy that can ignite giant emission. 

Also the case of Vela pulsar is very interesting in view of the models of giant emission. 
Kramer, Johnston & van Straten (2002) reported giant micro-pulses in this pulsar, which are 
different from normal micro-pulses and are probably closely related to classical giant pulses. 
They have typically large amplitudes, appear to be narrower than normal micro-pulses and 
exhibit a power law in their cumulative probability distribution. The first two properties are 
quite natural within our model (a sequence of unresolved nano-pulses?) and we discuss the 
last issue bellow. 

According to our scenario, both ordinary and giant pulses are related to the inner gap 
sparking activity and originate at relatively low altitudes, contrary to the suggestion that the 
latter arise in the outer gap region (Romani & Johnston 2001). Generally, the following two 
steps are involved: first a corresponding spark should be relatively intense (which corresponds 
to relatively high Lorentz factors), so the associated solitons are distinguished from the 
background radiation (the soliton radiation intensity depends very strongly on the Lorentz 
factor, see eq.[43] in GLM04). Then, in the second step, the intensity of the curvature 
radiation of distinguished soliton(s) can be kinematically boosted by means of equation (5). 
It is intuitive that the distribution of spark energy is quasi-gaussian, with a low energy cutoff 
corresponding to the pair formation threshold. Since the giant pulses should be associated 
with sparks corresponding to the high energy tail of this distribution, it is natural that they 
exhibit a power law in their cumulative probability distribution. 

This paper is supported in part by the Grant 2 P03D 008 19 of the Polish State Com- 
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